Does Singlet Fission Enhance the Performance of Organic Solar cells? 



Or 

CD 

c/3 



K. Aryan/pour, J. A. Munoz, and S. Mazumdar 
Department of Physics, University of Arizona Tucson, AZ 85721 
(Dated: September 21, 2012) 

Singlet fission, in which the optical spin-singlet exciton dissociates into two low energy triplet 
excitons, has been proposed as a viable approach to enhance the quantum efficiency of organic solar 
cells. We show that even when singlet fission is occurring in the donor molecule, the electronic 
structure at the donor-acceptor interface must satisfy specific requirements for the solar cell perfor- 
mance to be enhanced by this process. We focus on the pentacene-C6o solar cell and based on our 
calculations and available experimental data we conclude that there is not enough evidence that 
these requirements are met by the donor-acceptor interface here. We propose experiments that can 
determine whether the minimal requirement for enhanced performance driven by singlet fission is 
met in this and other solar cells. 



Introduction 

Multiple exciton generation (MEG), involving the 
generation of two or more low energy excitons from 
the absorption of a single high energy photon 1 -, has 
been suggested as an important means to overcome 
the Shockley-Qeisser upper limit of 33%2 for the 
quantum efficiency (QE) of inorganic solar cells. 
Within the MEG scenario, an optically generated 
high energy exciton undergoes conversion to several 
lower energy excitons, while obeying energy conser- 
vation. If each of the low energy excitons now un- 
dergoes dissociation into electron and hole carriers, 
the QE will exceed the limit for a single exciton. A 
related process, singlet fission (SF), has generated 
considerable excitement in the context of organic 
solar cells^r— . In organic 7r-conjugated materials, 
total spin is usually a good quantum number, and 
exchange interactions are large. The spin selection 
rule limits optical absorption to spin singlet states 
only, with most of the oscillator strength concen- 
trated in the lowest optical exciton in the quasi-one- 
dimensional material o 18 i 19 . It has long been recog- 
nized that in many 7r-conjugated systems, the ener- 
gies of the optical spin singlet state and the lowest 
triplet state satisfy the inequality Eg > 2Et ± , where 
Es (-Eti ) is the singlet (triplet) exciton energy. In- 
deed, in many such systems, the lowest singlet state 
is not the optical exciton, but an optically forbidden 
two-photon state that is an entangled state of two 
triplets (hereafter TT)22r— . In principle, such a sys- 
tem is a candidate for SF, whereby the singlet optical 
exciton undergoes fission into two triplet excitons. 
While SF in organic materials has been known for a 
long time^, recent excitement began with the obser- 
vation of relatively high power conversion efficiency 
of organic solar cells with pentacene (hereafter PEN) 
as the donor (D) molecule and Cgo as the accep- 
tor (A) 24 i 25 . Experimental demonstrations of SF 
in tetracene and pentacene crystals^— ' n i 16 ' 17 i 26 ~ — 



have led to the idea that the enhanced performance 
of PEN-C 60 solar cells is due to SF. 

With few exceptions 5 *^, research on SF has 
largely focused on the determination of the mech- 
anism of SF in the acenes2r-i~— i22i2£n2I . Whether 
or not SF can give enhanced performance, however, 
depends on the electronic structures of both D and 
A, and in particular, of the D-A interface. The goal 
of our work is different from the bulk of the existing 
theoretical work on SF— 10 ' 12 and is complementary 
to this research; it is to determine the conditions 
which need to be satisfied for SF-driven higher QE 
in organic solar cells. 

At the heart of organic photovoltaics lies pho- 
toinduced charge-transfer (PICT) between D and 
A at their interface. Optically excited D (A) do- 
nates an electron (hole) to its partner, generating a 
charge-transfer exciton D + A~. The lowest energy 
charge-transfer exciton (hereafter CTXo), depend- 
ing upon its binding energy, now undergoes recom- 
bination as well as charge-separation, and only the 
latter process is useful in photovoltaics. For SF to 
give enhanced photovoltaic performance, each of the 
two molecular triplets should now donate an elec- 
tron or a hole to its partner. One then sees right 
away that for SF-driven higher performance each 
of the following conditions have to be satisfied: (i) 
-E-CTXo < Et 1 , where i?cTX is the energy of the 
lowest charge-transfer exciton, (ii) the binding en- 
ergy of CTXo should not be prohibitively large, and 
(iii) the ground-state should continue to be neutral 
covalent and should not have undergone transition 
to an ionic state (such neutral-to-ionic transition, for 
example, occurs in crystalline mixed-stack charge- 
transfer solids-™.) Note that condition (i) requires a 
very low energy CTXo, given that Et 1 should sat- 
isfy Es > 2-Eti • Indeed, such low energy triplets are 



covalent in the valence bond languagi 



,21.22 



suggest- 



ing their unsuitability in charge-transfer processes. 
Conditions (ii) and (iii) have to be therefore satisfied 



in spite of the very low EctX • 

In the work presented here, we report explicit cal- 
culations of PICT on the PEN-Ceo system. We show 
from comparisons of detailed calculations and avail- 
able experimental information that whether or not 
SF in PEN is expected to give higher performance 
of the PEN-Ceo solar cell is not at all obvious. Even 
if SF-driven enhanced performance is occurring in 
this system, we show that PEN-Cgo is a marginal 
case where the above conditions are barely satisfied, 
and it cannot be assumed that all the molecules 9 
that are being investigated as candidates for SF will 
necessarily lead to higher QE for PICT. Elaborate 
evaluations of the above conditions for each D- A pair 
are essential, since for systems in which they are not 
satisfied, SF provides a competing channel for the 
decay of the photoexcitation and diminished perfor- 
mance. 

In the next section we present our theoretical 
model and discuss the computational approach we 
have taken to simulate the PEN-Ceo interface, and 
to determine the binding energy of CTX . Fol- 
lowing this, we present our results, for parameters 
appropriate for isolated PEN and Cgo molecules, as 
well as for the solid state heterostructure. Although 
our principal goal is to determine trends and to 
provide insight, we show that our calculations 
reproduce many of the known experimental results 
almost quantitatively. This is because our parame- 
terization of the PPP Hamiltonian has been done 
with considerable care^i. Our calculated absolute 
exciton energies as well as the exciton binding 
energies in single-walled carbon nanotubes 3 ^, and 
of one and two-photon states in polycyclic hydro- 
carbons that are molecular fragments of graphene 36 
have shown similar quantitative agreement between 
theory and experiment. We present the results for 
the energetics of the charge-transfer excitons for 
the PEN-Ceo interface, and its binding energy, for 
a geometry that should be maximally favorable for 
PICT. Based on these computational results, and 
additional work on the extent of charge-transfer 
in the ground-state presented in the Appendix, we 
present our conclusions in the final section. 

Theoretical Model and Methods. 

We have performed correlated-electron calcula- 
tions for an idealized quasi-one-dimensional het- 
erostructure with four PEN and three Ceo molecules 
forming a perfect one-dimensional stack (4PEN- 
3C6o hereafter, see Fig. [Ha)). We assume the PEN 
molecules to lie directly above or below one an- 
other with an eclipsed geometry. We assume sim- 
ilarly that one hexagonal face of each Cqq faces a 
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FIG. 1: (Color online) (a) Schematic of the PEN-Ceo 
"heterostructure" investigated in our work (see text); 
intermolecular separations are as indicated in the fig- 
ure, (b) The PPP HOMO and LUMO energies of PEN 
and Ceo within the HF approximation, and with our 
Coulomb interaction parameters; the Ceo HOMOs and 
LUMOs are five and threefold degenerate, respectively, 
(c) Calculated energies of the singlet optical exciton in 
PEN and the lowest PEN + C^" charge-transfer exciton, 
relative to the ground-state, using parameters appropri- 
ate for isolated PEN and Ceo, and identical screening pa- 
rameters for intra- and intermolecular Coulomb interac- 
tions. The energy of the triplet Ti is from experiment 28 . 
For SF to give enhanced solar cell efficiency, the lowest 
charge-transfer exciton should occur below Ti, as the 
curved arrow in the figure indicates. 



hexagonal face of the next Ceo molecule, with all 
intermolecular carbon atoms perfectly aligned. Fi- 
nally, a hexagonal face of the top Cgo molecule is as- 
sumed to be perfectly aligned with the central ben- 
zene nucleus of the proximate PEN molecule (see 
Fig. [TJa)). The minimum separations between two 
Cgo molecules as well as that between the closest 
PEN and Cgo are 0.35 nm in our calculations, while 
that between the PEN molecules are taken to be 0.40 
nm. These separations are representative of inter- 
molecular separations in organic molecular crystals 
(including, in particular, epitaxially grown films of 
Cgo on VSe^)- We are aware that the relative ori- 
entations between PEN and Ceo are quite different 
in the real system s 38 ' 39 , but have deliberately chosen 
this idealized conformation, since this will promote 
maximally stable CTX and most efficient charge sep- 
aration, because of the large intermolecular hoppings 
that result from the idealized geometry. While re- 
combination is also higher with this geometry^ -, this 
is not of concern here. 
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Our calculations are within the Pariser-Parr- 
Pople (PPP) Hamiltonia n 41 ' 42 for a two-component 
system 4 ^, 

H Hi n i ra -\- Hi n f er (1) 
The intramolecular component Hi ntra is given by, 

H intra = - t^(cl ia c^. a + H.C.) + 

Y £» n »,i,<? ■ ( 2 ) 

fj,— 5 — 7 

In the above c i a creates a 7r-electron of spin a 
on carbon atom i of molecule /i, with /i = 1 — 4 
corresponding to PEN, and fj, = 5 — 7 correspond- 
ing to C 60 , respectively, n^,^ = ^ CT c Mi i )CT is the 
number of electrons on atom i of molecule with 
spin tr and n^,j = J^ CT n^^^ is the total number of 
electrons on atom i of the molecule. U and Vij are 
the on-site and intramolecular intersite Coulomb in- 
teractions, respectively. The Vij are obtained from 
a modification of the Ohno parameterization 44 , 

Vij = [7/^1 + 0.6117^, (3) 

where By is the distance between carbon atoms i 
and j in Angstroms and k is an effective dielectric 
constant. Previous work has shown that U — 8 eV 
and k = 2 give excellent fits to absolute exciton en- 
ergies as well as exciton binding energies in the tt- 
conjugated polymer poly-paraphenylenevinylene 34 
and single- walled carbon nanotubes 3 -^. We have used 
standard nearest neighbor one-electron hopping in- 
tegrals t§ = 2.4 eV for phenyl C-C bonds in PEN 34 , 
and iy = 1.96 eV and 2.07 eV, respectively, for the 
bonds within the pentagons of Cqq and for those 
linking the pentagons. The smaller Ceo hopping in- 
tegrals reflect the curvature that reduces the overlap 
between neighboring p-orbitals 35 . We show below 
that excellent fits to various energy gaps of isolated 
PEN and Ceo are obtained with these parameters. 

The intermolecular component Hi nter of the 
Hamiltonian is written as, 

Hinter = ~ Y *tf ( C I,»>V^> + -ff-C.) + 

\ Y ^Ki-iXvri)' ( 4 ) 

We use the same functional form for VA as for 
the intramolecular Coulomb interaction with how- 
ever a dielectric constant k± that may be different 



from k (smaller n± implies greater intermolecular 
interaction). Intermolecular hopping integrals fol- 
low tjj — /3exp[(c — dij)/S], where the prefactor 
(3 = 0.2 eV, c is the minimum vertical distance 
between the molecules, dij is the distance between 
atom i belonging to molecule fi and atom j be- 
longing to molecule // =/= /i and the decay constant 
5 = 0.045 nm 4 ^. 

The last term in Eq. [2] contains a site energy 
e M that is nonzero only for the atoms of the Ceo 
molecules, and that is included to manipulate the 
energy offsets between the molecular orbitals (MOs) 
of PEN and Ceo, in order to vary the energy of 
CTXo. The energy of the charge-transfer exciton 
in an arbitrary D-A system, in the limit of zero in- 
termolecular hopping, is approximately given by 4 ^ 
I Ed — EAa — Eq , where I Ed is the ionization en- 
ergy of D, EAa is the electron affinity of A and 
Eq is the Coulomb stabilization energy due to prox- 
imate oppositely charged D and A in D + A~ . This 
concept in the past has been used mostly for molec- 
ular D and A with nondegenerate highest occupied 
molecular orbital (HOMO) and lowest unoccupied 
MO (LUMO). Very interestingly, we find that it is 
equally well applicable here, with fivefold degenerate 
HOMO and threefold degenerate LUMO in Ceo (see 
below). I Ed and EAa are both one-electron quan- 
tities, and their values in the gas phase and in the 
solid state can differ widely. We simulate the modi- 
fications of I Ed — EAa in the solid state empirically 
by the effective site energy e^, which modulates the 
mean-field Hartree-Fock (HF) energy difference be- 
tween the LUMO of the acceptor and HOMO of the 
donor, in accordance with Koopman's theorem. We 
shall refer to this energy difference as Alh in what 
follows (i.e., Alh = IEd - EAa). The Coulomb 
stabilization energy Eq on the other hand originates 
from the Vy" in Eq. |H 

Our calculations are mostly within the single 
configuration interaction (SCI) scheme, including CI 
between all one electron - one hole excitations from 
the HF ground-state, which we take to be a product 
function of the HF ground-states of the individual 
molecules. This approach enables us to determine 
ionicities of excited states quantitatively^,, and 
also the location of the excited electron and hole 
in a charge-transfer state. As we show below, this 
ability to precisely characterize all excited states 
allows us to determine the binding energy of CTXo 
for the model system we are considering. Finally, 
the SCI assumes the ground-state to be neutral and 
whether or not there is ground-state charge-transfer 
cannot be determined using this approach. We 
report separate full CI calculations on simpler 
model systems (see Appendix) to demonstrate 
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that conditions Eg > 2Et 1 and -Ectx < Eti can 
be simultaneously satisfied even with a neutral 
ground-state. 

Results 

We begin with computational results for the 
case where HOMO-LUMO offsets between PEN 
and Ceo are assumed to be the same as in the gas 
phase (e M = 0) and the screening parameter for the 
intermolecular Coulomb interactions is assumed to 
be the same as for the intramolecular interactions 
(k± = k). Following this we simulate solid state 
effects by performing calculations for nonzero e M 
and varying k±. Finally, we report calculations of 
the binding energy of CTXo as a function of e M and 

K±. 

PEN-Cf;n interface with "gas phase" Alh. In 
Figs. [Ub) and (c) we have given our results for 
isolated molecules. Fig. [IJb) shows the calculated 
HF HOMOs and LUMOs of PEN and C 60 . Our 
calculations are within the 7r-electron approxima- 
tion, and the absolute HOMO and LUMO energies 
are not meaningful. However, all energy differences, 
including the calculated Alh, are relevant. Fig.QJc) 
shows the main result for the PEN-Ceo interface, 
using the MO energies and wavefunctions corre- 
sponding to Fig. QJb) and k± = n = 2. Thus the 
results of Fig. QJc) correspond to the assumption 
that the molecular HF energies and wavefunctions 
are not perturbed at all at the interface in the solid 
state. Our calculations are for the singlet states 
only and the energy of the triplet exciton is taken 
from experiment^. 

There are several items of interest in the results 
shown in Figs. QJb) and (c). Our calculated energy 
of the optical exciton in PEN, 2.10 eV, is extremely 
close to the experimental energy of 2.07 eV in 
solution^. Our calculated HOMO-LUMO gap 
for C 6 o 3.37 eV is practically the same as the 
experimentally determined value of 3.36 eV— . 
The calculated energies for the two lowest allowed 
optical absorptions in the gas phase of Ceo, 3.1 
and 3.5 eV (not shown) are also very close to the 
experimental values 3.0 and 3.6 eV respectively for 
Ceo molecules dissolved in decalinii. All of these 
give confidence that our computational approach 
and the parameters used reproduce the behavior of 
the individual molecules almost quantitatively. Our 
calculated Ectx of 1.77 eV, taken together with 
the calculated Alh = 2.50 eV indicates that Ec is 
close to 0.7 eV, provided the assumption that the 
approximate expression i?cTX = Alh — Ec is valid 
holds here. This value of Eq is slightly larger than 



that calculated from electrostatic considerations 
for PEN + C^ at an intermolecular distance of 0.35 
nm (see Fig. 5 in Reference liol). Had we evaluated 
-EctXo from I En — EAa — Ec, using the known 
"bare" ionization energy of PEN 6.59 eV— and 
bare electron affinity of Ceo 2.68 eV— , along with 
the calculated Ec of Reference |40j, the calculated 
quantity would have been significantly larger than 
the 1.77 eV obtained by us. It is then reasonable to 
hold our calculated Ectx as a realistic lower limit 
for the energy of the charge-transfer exciton with 
gas phase parameters. As indicated in Fig. (He), 
enhanced QE due to SF is not expected here. 
The curved broken arrow indicates the extent to 
which the energy of CTXo needs to be lowered for 
SF to give higher QE. Smaller Alh and/or larger 
Ec than in the gas phase would be necessary for this. 

Simulation of solid state effects. Smaller Alh in the 
solid state is a consequence of the smaller ionization 
energy of PEN and the larger electron affinity of Cgo 
in the solid state. As mentioned in section II, we 
simulate solid state effects phenomenologically by 
varying e M . In Fig. [2] we show the calculated Ectx 
as a function of e M for several different /cj_ that enter 
into the calculations of Vjj- . The linear variation of 
Ectx with is surprising, as this indicates that 
the simple expression -Ectx — IEb — EAa — Ec 
continues to hold for a broad range of e M and k±, 
even for nondegenerate HOMO and LUMO in Cqq 
and with nonzero electron hoppings between PEN 
and Ceo- It is conceivable that this is unique to 
Cgo as the acceptor, since there can be few direct 
C-C intermolecular hoppings in this case. For each 
Kj_ there exists a critical site energy at which 
EctXo < Et x . 

The linear variations of Ectx in Fig. [5] indicate 
that Alh also varies linearly with e^. From the 
calculated e£ = 0.89 eV for k± = 2 (Fig. |2J) we 
find that the critical Alh at which Ectx (i < -^Ti 
is ~ 1.6 eV for the case where intra- and inter- 
molecular screening are comparable. We conclude 
that the lowest PEN + Cg charge-transfer exciton 
will be above the lowest triplet exciton unless at 
least one of the following two conditions are met (i) 
A LH < 1.6 eV (ii) |^7 C | > 0.7 eV. 

Binding energy of the lowest charge-transfer exciton 
Our definition of the binding energy of CTXo, Eh, 
is the usual one: it is the energy difference between 
CTXo and the lowest state in which the hole on 
the PEN molecules and the electron on the Cgo 
molecules are free. A relevant question to ask 
therefore is whether the very low energy of the 
CTX that satisfies Ectx < ^Ti implies also 
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FIG. 2: (Color online) -Bctx versus e M , the site energies 
on the carbon atoms of Ceo, for different k±. The dashed 
line parallel to the abscissa is the experimental Et 1 from 
Reference [2a. 



a large binding energy. If true, this would also 
impact solar cell performance negatively. We have 
calculated Et, for our model heterostructure as a 
function of e M to investigate this question. 

The calculation of E\, is nontrivial, since in spite 
of the multiple PEN and Ceo molecules retained in 
our calculations the overall system of Fig.QJa) is still 
discrete and there is no true continuum. Thus the 
binding energy cannot be determined from energy 
considerations. We determine the threshold of the 
continuum from wavefunction analysis. We identify 
a specific high energy CT exciton with widely sepa- 
rated and delocalized electron and hole as the lower 
threshold of the continuum, and calculate the energy 
difference between this state and CTX as the lower 
bound for ■ The wavefunction analysis is however 
complicated because in addition to states with com- 
plete charge-transfer there occur in the same energy 
range many other excited states including neutral 
Ceo optically dark states and states with incomplete 
charge-transfer 43 . We ignore these additional irrele- 
vant states and consider only excited states with at 
least 90% charge-transfer. 

In Fig. we show schematics of CTXo and 

the high energy charge-transfer exciton, hereafter 
CTX n , that we use in our evaluation of E^. 
The arrows indicate charge-transfer from the PEN 
molecules to the Cgo molecules identified in the fig- 
ure. The hole in CTXo is delocalized over two PEN 
molecules, with the electron on the nearest Ceo- 
The hole and the electron are both delocalized in 
CTX n which has relatively week contribution from 
configurations with the electron on the Cqq nearest 



to the PEN-Ceo boundary. The wavef unctions of 
these excitons are nearly independent of and k± . 
In Fig. Efb) we have plotted the energy difference 
between higher energy charge-transfer excitons and 
CTXo against the quantum number of the former, 
for kj_ = 2 and kj_ = 0.95, and for several different 
e M in each case. The terminal points for all curves 
correspond to CTX n . The consequence of smaller 
Alh (larger \e^\) is to shift the quantum number 
of charge-transfer states relative to the "irrelevant" 
neutral states, but has no bearing on their energies. 
The exciton binding energy E^ is independent of 
for fixed k± as is shown in Fig.[3tc) for several dif- 
ferent k±. This result has important positive impli- 
cation for the effect of SF on the solar cell efficiency: 
if low energy for CTXo is reached predominantly due 
to the reduction of the ionization energy of PEN 
and enhancement of the electron affinity of C^o in 
the PEN-Ceo heterostructure, the binding energy of 
CTXo is not affected and continues to be small. In- 
terestingly, our calculated binding energy for CTXo 
(0.46 eV) for k± — 2 is nearly identical to the es- 
timation by Verlaak et al. (0.438 eV) from electro- 
static considerations for the charge-transfer exciton 
in the PEN-Ceo heterostructure^. It is likely, how- 
ever, that in a real three-dimensional structure this 
is somewhat smalle r 50 ' 51 . 

The actual decrease in Alh due to solid state 
effects in the real system may not be sufficient to 
give Ectx < Et ± (see next section), in which case 
the lowering of -Ectx has to be driven by larger 
\Ec\- We have therefore calculated the critical k± 
required to obtain -Ectx < Et 1 for several different 
Alh > 1.6 eV, and for each critical k± we have de- 
termined E\j. Our calculated results for E\> against 
Alh ; now for different k± necessary to bring CTXo 
below Ti, is shown in the inset of Fig. Eld). We see 
that i?b increases dramatically and steeply as kj_ is 
reduced. 

To summarize: (i) PPP calculations within 
the SCI approximation indicate that the approx- 
imate expression IE-q - EAa - Eq for the lowest 
charge-transfer exciton is remarkably accurate for 
PEN-Ceo over a broad range of e M and k±; (ii) Low 
Ectx < -E-Ti in the solid state does not imply a 
high exciton binding energy of the CTXo exciton 
if the lowering of its energy was primarily due to 
the smaller Alh in the solid state; (hi) on the other 
hand, if the lowering of the CTXo energy is because 
of the larger intcrmolccular Coulomb interaction 
in the heterostructure, the binding energy will be 
substantially enhanced. Higher performance of 
the solar cell in the latter case is not expected, in 
spite of SF. Reviewing of the available experimental 
information, in particular of Alh, thus becomes 
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FIG. 3: (Color online) (a) Schematic of the lowest energy charge-transfer exciton CTXo. The arrows denote the lengths 
of the charge-transfer "bonds" . The wavefunction of CTXo is predominantly a superposition of the configuration with 
the hole on the PEN molecule and the electron on the Ceo molecule at the PEN-C60 interface and the configuration 
with the hole on the neighboring PEN molecule, (b) Schematic of the high energy charge-transfer exciton CTX n 
with delocalized hole and electron. The arrows have the same meaning as in (a), (c) Energy differences between all 
higher energy charge-transfer excitons and CTXo, versus the quantum numbers of the higher excited states. The 
terminal point on each plot is CTX„ of Fig. EJb). a — k±/k\\ is a normalized measure of the inter- over intraunit 
screening, (d) Eh versus e M for three different a. The inset shows Eh against Alh, each with a values required to 
satisfy -Ectx < Et x . 



essential for determining whether or not SF is 
beneficial. 

Discussions 

As mentioned in the Introduction, the bulk of 
the existing literature on SF is on single-component 
PEN, with the focus on understanding the mecha- 
nism of SF. To the best of our knowledge only three 
different groups have investigated the PEN-Cgo het- 
erostructure in the context of SF— i 14 i 16 , and have 
concluded that SF indeed enhances charge genera- 
tion. We reexamine aspects of these investigations 
carefully below, while emphasizing that our goal is 
not to claim that their conclusions are incorrect, 
but rather, to point out that further experimental 
and theoretical work is necessary to establish be- 
yond doubt that SF is indeed behind the relatively 
large QE of PEN-C 60 solar cells^s. 

Among the works of interest, only the work by Lee 
et al. has made independent attempt to evaluate 
the QE of charge generation: the authors claimed 
45% enhancement of the external QE for photocur- 
rent generation in a PEN-Ceo multilayer photode- 
tector, and ascribe this to enhanced internal QE due 
to SF in PEN 6 . This work is not based on the so- 
lar cell configuration, and it is not clear whether 
it is possible to predict solar cell performance from 
it. Furthermore, the estimation of the photocurrent 
efficiency was dependent upon various assumptions 



and appears to be model dependent. Weak magnetic 
field dependence of the photocurrent upon laser ex- 
citation close to PEN absorption wavelength was ob- 
served, and was explained as due to reduction in the 
singlet character of the TT state in the presence of 
magnetic field. The TT singlet state is, however, 
a quantum mechanical eigenstate of the full many- 
body Hamiltonian and it is not clear why the nature 
of this wavefunction should change in the presence 
of a weak magnetic field. We conclude that while 
this work is of significant interest, additional inde- 
pendent measurements would be necessary to sub- 
stantiate the principal claim of this work. 

Rao and Wilson et a/. 16 ' 17 , and Chan et alr^ 
actually have proposed very different mechanisms 
of electron transfer from the excited PEN state. 
Rao and Wilson et al. 's interpretation of their ul- 
trafast spectroscopic measurements is the standard 
one, and the authors do believe that CTX is be- 
low the triple t 16 ' 17 . This latter conclusion is based 
upon an estimate of 4.5 eV for the EA of Cgo in the 
solid stated, which is larger by 1 eV compared to 
other estimates 37 . Such a large EA for Cgo, taken 
together with the solid state ionization energy of 
PEN (5.1 eV)££ give a sufficiently small A LH that 
could actually lead to possible neutral-ionic transi- 
tion in the ground state (see Appendix). Delayed 
charge generation is ascribed to longer triplet dif- 
fusion time to the heterojunction interface. Note, 
however that longer diffusion time is also expected 
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if the excited state in PEN is a superposition of TT 
and a PEN excime r 53 i 54 . Such a superposition in 
fact is suggested from theoretical work— Standard 
one-electron transfer from the excimer component of 
the wavefunction is expected in this case, with the 
high QE (but less than one!) due to the long lifetime 
of the excimer and the ease of hole derealization in 
PEN. 

Chan et a/— have contradicted the scenario of 
single electron transfers from individual triplets and 
have proposed multiple electron transfers from either 
a quantum mechanical superposition of the PEN sin- 
glet exciton and the TT state (the authors refer to 
this last state as a multiexciton state and denote it 
as ME), or from a related multiexciton state which 
they label as ME'. The possibility that the ME 
or ME' state donates only one electron to one Cqq 
molecule, is discarded by Chan et al. based entirely 
on the earlier work by Lee et al£. Note that ME 
is also expected to have an excimer component^, 
which would make single-electron transfer from the 
excimer component of ME at least a likely chan- 
nel. Assuming that the TT component of ME or 
ME' involves more than one PEN molecule^ the 
proposed multi-electron transfer mechanism also re- 
quires i?cTX < £ti, since the multiexciton state 
needs to generate two charge-transfer excitons. It is 
this last requirement that we now examine below. 

Based on our computational results, we con- 
clude that the upper limit for Alh that can give 
SF-enhanced performance is 1.6 eV. This is justi- 
fied from the demonstration that Et, rises sharply 
for larger Alh- Note also that Eq calculated 
by others^ is smaller than found by us, which 
would necessitate even smaller Alh- Coinciden- 
tally, the ionization energy of PEN in the solid state 
(5.1 eV)£^, taken together with the electron affin- 
ity of Ceo as measured by Schwedhelm et al. (3.5 
eV)2£ gives Alh exactly equal to this critical value. 
It has, however, been demonstrated that the en- 
ergy mismatch between levels in heterostructures 
cannot be obtained from studies of the individual 
semiconductors— We are aware of only two refer- 
ences that have attempted to determine Alh based 
on experiments on the PEN-Ceo heterojunctio n 56 ! 57 
(as opposed to estimates based on experiments on 
isolated systems.) Kang et aZ— find Alh of 1.56 eV 
and 1.50 eV, respectively, for Coo deposited on PEN 
(with gold as the substrate), and PEN deposited on 
Ceo (also with gold as the substrate). Thus in both 
cases the measured Alh is considerably smaller than 
what would be obtained from EA of 4.5 eV for Cgcr^ 
assumed in Reference [lfl and would barely satisfy 
the necessary condition for enhanced QE. Signifi- 
cantly different Alh are found by Salzmann et aZ.j5i 



who find this quantity to depend strongly on the pro- 
cessing technique used to generate the heterostruc- 
ture. For layered structures of Cgo on PEN precov- 
ering PEDOT:PSS the authors find A LH ~ 1.15 eV 
(see Fig. 3(a) of Reference H3), which would cer- 
tainly make enhanced QE feasible. However, for 
codeposited films of PEN and Cgo the authors de- 
termine Alh ~ 1-75 eV (see Fig. 3(b) of Refer- 
ence H3), which would make enhanced QE unlikely. 
Even more importantly, in both references only the 
HOMO energies of PEN and Ceo are determined di- 
rectly from experiments. The LUMO energy of Ceo 
is estimated from the transport gap— The latter ap- 
proach makes the estimate of Alh questionable, as 
recent work has shown that the transport in molec- 
ular solids involves mteraiolecular charge-transfer 
states that can occur below the molecular LUMO 
level— In the present case, while References l56l and 
l57l assume the HOMO-LUMO gap of solid C 60 to be 
2.6 eV, the same quantity is estimated to be 3.36 
eV for epitaxial layers of Ceo on layered VSe2 sin- 
gle crystals by Schwedhelm et a/—, and even larger 
figures had been given by earlier investigators^— . 
The HOMO-LUMO gap of Schwedhelm et al, taken 
together with the HOMO offsets between PEN and 
Ceo determined experimentally in References I56H571 
would put Alh outside the region that could give 
enhanced QE. 

In summary, we have presented a careful the- 
oretical analysis of the PEN-Cgo interface within 
the correlated 7r-electron PPP model that is able 
to reproduce many of the experimental features of 
PEN and Ceo quantitatively. We simulate solid 
state effects by independently varying parameters 
that change the energy separation Alh between 
the LUMO of C 60 and the HOMO of PEN, and the 
many-body Coulomb interaction that contributes to 
the binding energy of the PEN + Cg~ charge-transfer 
exciton. Based on these calculations we conclude 
that while it cannot be ruled out that SF is behind 
the high QE of the PEN-Cgo solar cell, neither is 
there unqualified support for this scenario from 
currently available experimental information. A 
variety of experiments have recently detected the 
charge-transfer exciton below the optical gap of 
the donor polymer in organic heteroj unctions^— 
Interestingly, in all such cases the charge-transfer 
exciton has been found to occur at 1.3-1.6 eV, viz., 
significantly above the 0.9 eV where the PEN + Cg~ 
charge-transfer exciton needs to occur for enhanced 
QE. We propose similar measurements for the 
detection of CTXo in PEN-Cgo- In particular, 
direct photoexcitation of CTXo has been possible 
in spite of its low oscillator strength— We propose 
that experiments be performed on the PEN-Ceo 
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FIG. 4: (Color online) (a) and (b): Interacting 
pairs of ethylene-like and butadiene-like molecules, with 
Coulomb parameters chosen such that the lowest molec- 
ular triplet is nearly at half the energy of the optical 
singlet. MO offsets created artificially generate D-A sys- 
tems in both cases, (c) and (d): i?cTX and the ground- 
state ionicity pgs versus the MO offset, for the systems 
in (a) and (b), respectively. 



hetero junctions to detect the charge-transfer exci- 
ton directly. 
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Appendix: Higher Efficiency versus Ground- 
State Charge- Transfer 



The SCI using HF basis assumes closed shell MOs 
as the ground-state. It is conceivable that as Alh 
continues to decrease, this assumption breaks down 
and the ground-state of the D-A system becomes 
ionic, as has been observed in mixed-stack crys- 
talline charge-transfer solids-^. We have considered 
this possibility because a triplet state that is at 
half the energy of the singlet exciton is covalent in 
the valence-bond languag o 21 i 22 , and it is conceivable 
that before the ionic CTXo can be below the co- 
valent Ti the ground-state itself becomes predomi- 
nantly ionic, which would have negative consequence 
for solar cell efficiency. Whether or not this happens 
cannot be checked within the HF approximation. To 
test whether the condition -EctXo < E^ 1 can be sat- 
isfied without the ground-state becoming ionic we 
have done exact diagonalizations of our Hamilto- 
nian for coupled hypothetical small molecules (see 
Figs. IDJa) and (b)). Our systems consist of pairs 
of molecules that resemble ethylene and butadiene, 
with the difference that we choose Coulomb pa- 
rameters such that the triplet exciton is nearly at 
half the energy of the singlet optical exciton, while 
keeping all intramolecular and intermolecular elec- 
tron hoppings realistic. This is achieved by tak- 
ing U — 5.56 eV and k±_ = n = 1.7 (which is 
not too far from realistic parameters). We include 
the site energy term in Eq. [2] for one of the two 
molecules in each case, thus breaking the symmetry 
and creating a D-A system. We keep varying Alh 
slowly and monitor both £"ctx an d the ground- 
state ionicity (the relative weight of the ionic con- 
figuration D + A~ in the wavefunction) pes- Our 
results for the two cases are shown in Figs. @Ic) 
and (d). In both cases the ground-state continues 
to be covalent where -Ectx < -Eti is reached, and 
the neutral-ionic transition occurs at a much smaller 
Alh- We have confirmed that this is true for other 
k±_ also. Essentially, as long as the intermolecular 
hopping integral has a weak role, as is necessarily 
true in real D-A systems with Cgo as the acceptor, 
< Et, can be reached even with covalent 



E, 



CTXo 



ground- state. 
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